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Abstract 

Double-well optical lattice in an insulating state is considered. The influence of 
atomic vibrations and mesoscopic disorder on the properties of the lattice are stud- 
ied. Vibrations lead to the renormalization of atomic interactions. The occurrence of 
mesoscopic disorder results in the appearance of first-order phase transitions between 
the states with different levels of atomic imbalance. The existence of a nonuniform ex- 
ternal potential, such as trapping potential, essentially changes the lattice properties, 
suppressing the disorder fraction and rising the transition temperature. 
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1 Introduction 



Optical lattices provide a convenient tool for modelling various many-body periodic struc- 
tures [1-4]. The attractiveness of optical lattices is caused by the possibility of changing 
their properties in a required way. Thus, the lattice period, the depth of potential wells at 
lattice cites, and the lattice geometry, all these characteristics can be rather easily varied in 
experiment by applying the necessary laser beams. Also, the strength of atomic interactions 
can be regulated in a wide range by the Feshbach resonance techniques [5,6]. Respectively, 
it is feasible to create many artificial periodic structures with the properties necessary for 
different applications. 

A special class of optical lattices is represented by the double-well lattices, where each 
lattice cite is formed by a double-well potential [7-16]. Such lattices, similarly to single-well 
lattices, can be in superfluid or insulating state. In addition, they can exhibit other states 
that cannot arise in the single- well counterparts. For example, they can exhibit the ordered 
states with nonzero mean atomic imbalance and with phase transitions between the ordered 
state and a disordered state with zero average atomic imbalance [7-16]. Another unusual 
property of the double- well lattices is the possibility of forming mixed states, where meso- 
scopic regions of the ordered phase are intermixed with mesoscopic regions of the disordered 
phase [17]. In a bulk system, or in a sufficiently large trap, these regions are randomly 
distributed in space. While, for atoms in a tight trap, the region distribution is governed 
by the trap geometry. The latter would be similar to the radial distribution of superfiuid 
fraction in a trap with a periodic crystalline-type atomic configuration [18]. 

The mixed state of a double-well optical lattice, with mesoscopic regions of different 
phases, exhibits a variety of interesting properties, distinguishing such a state from the purely 
ordered or completely disordered states. In our previous publication [17], we considered a 
rigid lattice with no external fields acting on the system. The aim of the present paper is 
twofold: to generalize the consideration by including external fields, and also to analyze the 
role of atomic vibrations in a double-well optical lattice. In Sec. 2, we derive a model of 
an insulating optical lattice in order to explain the physical meaning of all its terms and 
for explicitly demonstrating the origin of the term representing the influence of external 
potentials. In Sec. 3, atomic vibrations are introduced and their role in the renormalization 
of atomic interactions is analyzed. In Sec. 4, we study a double-well optical lattice with 
mesoscopic disorder in the presence of an external field. Section 5 contains conclusions and 
discussion. 

Throughout the paper, we use the system of units with h — 1 and ks — 1- 



Deriving the model of a double-well optical lattice, we start with the standard Hamiltonian 



where ■0(r) is a field operator of atoms, $(r) is a pair interaction potential, and in the 
single- atom Hamiltonian 



2 Insulating Double- Well Optical Lattice 




(1) 




(2) 
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the external potential is a sum 



U{r) = 14(r) + U,,t{r) , 

where the first term is the lattice potential, periodic with the lattice spacing a, while the 
second term is not periodic, representing an additional external field, e.g., corresponding to 
a trapping potential. Atoms can be either bosons or fermions. 

Keeping in mind an insulating lattice, the field operator can be expanded over an or- 
thonormal basis of locahzed orbitals, 

H^) = J2cnjMr-rj) , (3) 

nj 

where ?7, is a quantum multi-index and rj is a vector of atomic location. With this expansion, 
Hamiltonian (1) transforms to the expression 

H — \^ \^ pmn J „ . I 1 ^nin2n3n4 t r ■ r ■ (A\ 

ij mn {j} {n} 

in which 

= / ^-(r - ri)H^{r)Mr - r,) dr (5) 

and are the matrix elements of the interaction potential. As localized orbitals, one 

can take maximally localized Wannier functions designed for double- well lattices [19]. 

Considering the lattice with unity filling factor, it is necessary to impose the no-double- 
occupancy constraint 

^ y ^lij^nj — 1 5 CfnjCnj — . (6) 

n 

And in the case of an insulating lattice, the no-hopping condition 

^mi^nj — ^ij^lnj^nj (7) 

is valid. Then Hamiltonian (4) becomes 

H — \^ \^ prnn J „ . , 1 T/nin2n3n4 t J r- -r ■ (R\ 

~ / ^ ii '-mj^nj ^ 2 / J / J ^ ij ^n\i^n2j^'^-i3^'^Ai i \° ) 

j mn i^j {n} 



in which 



^nin2nzn4 _ ^nin2n3n4 i ^nin2n3n4 
ij ~ ijji ijij 



the upper sign corresponding to Bose and the lower, to Fermi statistics. In the sums over 
atomic locations, one has i ^ j- To simplify the notation, we shall omit in what follows this 
inequality, setting instead the condition 

y-nin2n3"4 — g 

Each lattice site is formed by a double-well potential. The hopping between difi'erent 
lattice sites is absent in an insulating lattice. But the tunneling between the wells of a 
double-well potential in the same site is, generally, present. Therefore, to take into account 
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the tunneling, it is necessary to consider minimally two energy levels. In what follows, we 
take two such lowest levels with n = 1,2. The necessity of taking at least two quantum 
states distinguishes the case of the double- well lattice from that of a lattice with single-well 
sites. The usual situation is when the ground-state function is symmetric with respect to 
spatial inversion, while the next-level state is antisymmetric, so that 

V^i(-r) = ipi{r) , '02(-r) = -V'2(r) . 

Because of this property, the off-diagonal matrix elements, such as V^]^^^ or V^'^^^, are exactly 
zero. Even if the above symmetry property does not hold, such off-diagonal matrix elements 
are usually much smaller than the diagonal ones of the type V^^^^^, V^'^'^^, V^^"^^ ^ and V^]-^^^. 
The nonzero matrix elements enter the following formulas through the combinations 

The latter are symmetric with respect to the change of the indices: 
Introducing the notations 

U^^ = I C(r - r.)[/(r)^„(r - r,) dr , (10) 
and invoking the symmetry properties, we have 

X""" = 5 X"" fill 

which gives 

EJ]^ = 5mnKJP + . (12) 

The matrix elements of the kinetic energy can be represented as 



Also, we define 



where 



The quantity 



E«^lj:{ES + e!!) = j:i; + u„N, (14) 



n ■ jjy ^2m 

3 3 
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Hj^-E^^-E]j^-Ujf-US (16) 
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plays the role of an external field acting on the atom in a j-th lattice site. And 

i 

has the meaning of the tunneling frequency characterizing atomic tunneling between the 
wells of a double well. Depending on the shape of the double-well potential, this frequency 
can be varied in a wide range [20] . 



Then we use the operator transformation 



c\.c2j = - iS] , cljcij = + iS] . (18) 



The operators 



1 / + + \ ^„ i 

2 



Sj=l {4jC2j + i-cy) (19) 

satisfy spin algebra, independently from the statistics of atoms, whether Bose or Fermi. 
Because of their spin algebra, these operators are called pseudospin operators. With these 
notations, Hamiltonian (8) reduces to the form 



(20) 



According to their properties, the operator describes atomic tunneling between the 
wells of a double-well potential in a j-th lattice site. The operator S^j corresponds to the 
Josephson current between the wells. And the operator S'^ characterizes atomic imbalance 
of these wells [14-17]. 

3 Vibrating Double- Well Optical Lattice 

Atoms, forming the lattice, are, of course, not immovable, but can vibrate around their 
lattice sites. In this section, we study the role of such vibrations. 
Each atomic vector Vj can be represented as the sum 

^3 = + (21) 

of the mean atomic location 

a, ^ (r,) (22) 
and an atomic deviation from the mean location, such that 

(u,) = . (23) 

The latter condition follows directly from definitions (21) and (22). 
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The interaction matrix elements (9) are functions of the atomic vectors, so that, for 
instance, 

= A{r,j) {Tij = ri- rj) , (24) 

and, similarly, all other matrix elements (9). Since we are considering an insulating lattice, 
the atomic deviations are supposed to be small, because of which the interaction matrix 
elements can be expanded in powers of these deviations, limiting the expansions by the 
second-order terms, e.g., as 



a/9 



U, 



Aij — A{a4j) + ^ ^ijUij - - ^ Afj 
where the relative deviation is 
and the notations are used: 

{aij = aj — aj) , 
d^A{sL, 



(25) 



(26) 



Aa — 



dA{a,,) 
da- 



da?dd 



In order to guarantee the validity of condition (23), one has to introduce the grand 
Hamiltonian 

i/ = ^-^A,-u,-, (27) 



with the Lagrange multipliers cancelling in the Hamiltonian the terms linear in u". Then, 
expanding the matrix elements in powers of the relative deviations (26), we come to the grand 
Hamiltonian 



j ^ ^ ij \ a/3 



+ 



a/3 



ox QX 



/fa. 



V, 



a/3 a /3 
ij '^ij'^ij 



a/3 



13 '^ij'^ij^ 



QZ QZ 



(28) 



To make the problem treatable, one needs to decouple the vibrational and atomic degrees 
of freedom. We use the following decouphng: 



(29) 



The existence of atomic vibrations leads to the renormalization of atomic interactions. 
In turn, atomic interactions renormalize the phonon matrix 

a/3 ^a/3 , Q ryoili I QX qx\ n 7-a/3 / 



_ A^f + 2B^'{S^S^) - 2It[{StS^) . 
The renormalized atomic interactions are 



(30) 



a/3 



j ) 



a/3 

Also, we introduce the notation 



(31) 



ij aP 



in which 



A 



(32) 



(33) 



^3 



As far as v4 is a constant, we have 



dapa'^j 



ij i j 



In the phonon matrix (30), the main term is Aff . Therefore one has 



(34) 



a 13 



1-3 



ap a p 
ij ^i ^j 



(35) 



Thus, Hamiltonian (28) can be represented as the sum 

H^Eo + Hph + Hat, (36) 
in which the first term is given by Eq. (32), the second term is the phonon Hamiltonian 



j ij a/3 

and the third term is the atomic Hamiltonian 

Hat^-J2 i^oS] + HjS]) + J2 {BijSfSJ - hiStsf) 

j ij 

The phonon spectrum uiks is defined by the eigenvalue problem 



(37) 



(38) 



jP 



' ^ ^ks~ ^ks^ks ) 



(39) 



where Bks is a polarization vector, k is quasi- momentum, and s, polarization index. The 
phonon operators bks are introduced through the transformation 



Pi = - -1= E V'miOks eks (bks - &Lfe,) 

ks 

% = ^= E (hs + bKg) 



t 1 „ika, 



(40) 



Using this, the phonon Hamiltonian (37) becomes diagonal, 

Hph^J2'^''4^lshs + \) ■ (41) 



ks 

This allows us to find the deviation correlation function 

/3 



ks 

Since in the efi^ective atomic interactions (31), we have i ^ j, they can be written as 

S,, = S(a,,) - J2 Btfiu-u^) , h = Aa.,) - ) • (43) 



(44) 



And the nonoperator term (32) becomes 

^ ' ij a/3 

The equation for the phonon spectrum (39) can be rewritten as 



^ks 



J Oil: 



For illustration, let us consider a lattice with cubic symmetry. Then we can define the 
effective spectrum through the average 



^1 = \Y.^Is. (46) 



where d is space dimensionality. This yields 



' = - - ^ A,e*-^*^ , (47) 



" m 



with the dynamical matrix 



a=l 

In the long- wave limit, using the property 
we have 



D^J^-li2^?r■ (48) 



^^^^E^«^(k-ai.)' (k^O). (49) 



2m 



For a cubic lattice, the long-wave phonon spectrum reduces to 

Wfe ~ c/c (A; ^ 0) , 

with the sound velocity 

c = 

and the isotropic dynamical matrix 



Da2 
2m 



(50) 



(51) 



(52) 



Here a is the lattice parameter of the cubic lattice. 

It is important to emphasize that phonons can arise only in the presence of long-range 
atomic interactions. For example, in the case of local interactions [4] , we have 

C/exp (-g^aj) , 

where 

and U and Qq are positive real parameters. This leads to 

Then the dynamical matrix (48) is negative: 

Dij ^ -2ql<^>ij ($ij > 0) , 

which makes it impossible to define a real phonon frequency. 
If is not Gaussian, but simply exponential, as 



U exp(-qo • aij) {U > 0) , 



then 



and the dynamical matrix is positive: 



Dij ^ 2 ^o^^j ^^^^ > °) ' 



making well defined a real phonon frequency. 
When the interactions are of power law, as 

U 



^ij ~ — (n > 0) , 



then 



'^3 r.'i 



a: 



The latter, for nearest neighbors, gives 



^ fn + 2 \ 

Real phonon spectrum exists for a positive dynamical matrix Dij. This means that the 
power law requires the condition 

n> d-2 . 

If the lattice is sufficiently large, one can resort to the replacement of the summation 
over momenta by the integration over the Brillouin zone B, so that 

1 ^ r dk 

V JsWY' 

Then, from Eq. (43), we have 

/ a I3\ f 1 

This defines the mean square deviation 

d 

^ J2 ) = {n>l)d . (54) 

a=l 

In this way, the effective atomic interactions (43), for a cubic lattice, become 

B^J = B{^,) - I 5^ , 4 = /(a,,) - I 5^ I- , (55) 

a a 

which shows how atomic vibrations change the interaction of atoms. 

4 Double- Well Lattice with Mesoscopic Disorder 

Atomic vibrations can be called microscopic fluctuations, since they are related to the fluc- 
tuations of separate atoms around their lattice sites. As is shown in the previous section, 
such fluctuations renormalize the strength of atomic interactions. 

There exists another type of fluctuations, whose occurrence can lead to a much more 
noticeable change of system properties. These fluctuations are termed mesoscopic, because 
each fluctuation involves many atoms that move coherently. The fluctuations arc called 
mesoscopic, since their typical linear size is much larger than the mean interatomic distance, 
but smaller than the system size. These fluctuations represent strong fluctuations of order 
parameters, as a result of which the system can be treated as a fluctuating mixture of different 
thermodynamic phases, say, an ordered and a disordered phase. Then the appearance of 
such mesoscopic fluctuations of competing phases can be understood as the occurrence of 
mesoscopic disorder. This type of fluctuations can lead to an essential change of system 
properties, as has been demonstrated for several kinds of matter of different physical nature, 
e.g., for ferromagnets [21,22], antiferromagnets [23,24], ferroelectrics [25], crystals [26], and 



dk 

coth — - - — -7 
\2TJ {2ttY 



(53) 
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other materials reviewed in Ref . [27] . A somewhat close behavior happens for the clustering 
nuclear matter [28,29]. 

In our previous paper [17], we have considered the influence of such mcsoscopic disorder 
on double-well optical lattices, when there are no external fields. In the present section, we 
generahze the consideration to the case of a nonzero external field. The latter can be caused, 
e.g., by the presence of a trapping potential. We study a fiuctuating mixture of two phases, 
an ordered phase, where there exists an average large atomic imbalance between the wells 
of double-well potentials, and a disordered phase, where such an imbalance is small. In the 
absence of external fields, the atomic imbalance of the disordered phase is exactly zero. 

After averaging over various configurations, produced by mesoscopic fiuctuations, as de- 
scribed in Refs. [17,27], we have, instead of Eq. (20), the effective Hamiltonian 

H = H,^H2, (56) 

in which each term 



(57) 



corresponds to the related phase labelled by the index u = 1,2. Each Hamiltonian part (57) 
acts on a weighted Hilbert space Hu- And the total Hamiltonian (56) is defined on the fiber 
space 

n^ni(^n2. (ss) 

The geometric weights Wi, of the phases are defined as the minimizers of the thermodynamic 
potential 

F = -TlnTre"'^^ , (59) 

under the normalization condition 

Wi + W2 — 1 , < < 1 . 

The quantities of interest are the average values of the tunneling intensity 



^ E (-^l)^ ' (60) 



N 



Josephson current 

and the mean atomic imbalance 



IE(^J)- (61) 



N 



Sty = 



^ E i^^)^ ■ (62) 



N 



Here the index u in the statistical averages implies the averaging with respect to Hamiltonian 
(57). The atomic imbalance (62) plays the role of an order parameter distinguishing an 
ordered and a disordered phases. This definition implies the constraint 

si > S2 . (63) 
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When external fields are switched off, the order parameter of the disordered phase has to be 
zero, which assumes the limiting condition 

lim S2 = . (64) 

To calculate the above quantities (60), (61), and (62), we need to invoke some approxima- 
tion for the Hamiltonian terms containing pairs of pscudospins. For this purpose, we resort 
to the Kirkwood approximation [30], according to which for a pair of operators, acting on 
the space Tii,, one has the decoupling 

S^S^ = g^j {{Sr).Sf + Snsf). - {Sr)ASf).) (i ^ 2) , (65) 

in which g\- is a pair correlation function and the index v at the angle brackets means the 
statistical averaging with the related Hamiltonian (57). This decoupling differs from the 
mean-field approximation by the presence of the pair correlation function. 

In the formulas below, we shall employ the notations for the average interaction intensities 

^ ^ E ^'i^^^i ' ^ ^ E ^F^.- • (66) 

We assume that the trap is sufficiently wide, so that the quantities 

n- = n Hj = Ho (67) 
can be treated as independent from the lattice indices. Also, we introduce the notations 



E^^w^Eo + wi-\A- -f< + fstj (68) 
and 

= w^Q - wlB^Xy , = WyHQ + wll^Sy . (69) 

Then Hamiltonian (57) reduces to the form 

3 

For the thermodynamic potential (59), we find 

F^F, + F2, F,^K- NTln f2cosh ^] , (71) 



2T 



where 



j, = ^ni + i2. (72) 

The tunneling intensity (60) becomes 

X, = % tanh ( ^] . (73) 



J. \2T 
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The Josephson current (61) in an equilibrium system is zero, 

1/. = . (74) 

And the atomic imbalance (62) reads as 

~ ~T' tanh ( -^^^ J • (75) 



In order to estimate the relation between the different interaction terms (9), we can take 
into account that the values V^^^^, K^^^^ and Ki^^^ are close to each other, so that we can 

''J '■J ''J 

set 

T/llll _ T a2222 _ T/1221 _ -ir 
^ij — ^ij — ^ij — ^ij ■ 

As a result. 

We shall use this simplification in what follows. 

To accomplish numerical calculations, we need to define the pair correlation function. 
Here, for simplicity, we set 

9ij - 1 , (76) 

which reduces the Kirkwood approximation to the mean-field form. For the convenience of 
numerical calculations, we define the dimensionless thermodynamic potentials 

in which 

It is also convenient to introduce the dimensionless parameters 

u=j, oj=j, h= — . (79) 
Then for quantities (69) and (72), we have 

^ = w,,uj , Y = ^^/^^ ' Y = + '^l^'^ ' (80) 

where 



= ^/{h + w,s,)^ + uj^ . (81) 
Measuring temperature in units of /, for the tunneling intensity (73), we get 

while the atomic imbalance (75) becomes 

= ^^±^ tanh ( . (83) 
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In the dimensionless notation, the thermodynamic potential (77) reads as 



2 cosh 



2T 



(84) 



(85) 



with 

K 

NI ' 4 

Recall that, in addition to the above equations, we have to define the geometric weights 
of the ordered and disordered phases as the minimizers of the thermodynamic potential G. 
To take explicitly into account the normalization condition, we introduce the notation 

w = wi , W2 — 1 — w . (86) 

Minimizing the thermodynamic potential (84) with respect to w, we find the equation 



w — 



2u + uj{xi - X2) + h{si - S2) - si 
4u — sj — si 



(87) 



defining the weight of the ordered phase. The weight of the disordered phase, according to 
relation (86), is 1 — w. 

When the gap between the energy levels Ejj and defined in Eq. (5), is small, as 
compared to the interaction intensity (78), then, in view of definition (17), one has a; <^ 1. 
In that case, we have 



Xj. = (cu ^ 0) 



Sy = tanh 

Therefore, the atomic imbalance for the ordered phase becomes 

' wh + uP'Si 



2T 



Si = tanh 



2T 



(88) 



The atomic imbalance for the disordered phase is 

^1 — w)h + (1 — wys2 



S2 — tanh 



2T 



(89) 



And the equation (87) for the geometric weight of the ordered phase takes the form 



w — 



2u + h{si — S2) — si 
4u — si — si 



(90) 



These are the main equations we need to solve, keeping in mind constraints (63) and (64). 
In addition, we have to compare the thermodynamic potential of the mixed system, (84), 
with that of a pure ordered system, when w — 1 and there exists just one order parameter 
Si given by the equation 

Si = tanh I 
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Note that there always exists a solution with w = 0.5 representing a degenerate case, when 
the competing phases are asymptotically equivalent. Among all possible solutions, we have 
to select the most stable one, which minimizes the system thermodynamic potential and 
satisfies all imposed constraints, such as Eqs. (63) and (64) and the condition < w < 1. The 
solution w — makes the value Si undefined, hence constraint (63) inapplicable. Therefore 
this case is unphysical and has to be excluded. 

The results of numerical calculations are presented in Figs. 1 and 2. The parameters u 
and h have been varied in a wide range. When u < 0, the most stable solution corresponds 
to the pure system, with no mesoscopic disorder and = 1 for all temperatures. Positive 
values of u allow for the appearance of mesoscopic disorder. According to definition (79), the 
positive parameter u can be interpreted as the ratio of the disordering interaction strength 
A to the ordering interaction strength I. Thence this ratio u can be called the disordering 
parameter. It turns out that for u > 0.5 the solutions are not continuous. This could mean 
that such a strong disorder parameter prohibits the states with mesoscopic fluctuations, at 
least in the frame of the present model. To consider the higher values of u requires to invoke 
a more sophisticated approximation than the mean-field one. Varying u in the interval 
< u < 0.5 shows that at T = 0, we have Si = 1 and the most stable is the state with 
w = 1. With rising temperature, the state with w = 1 and S2 = persists till the transition 
temperature Tq, where there occurs a first-order phase transition to the state with w = 0.5. 
The transition temperatures, as functions of u and h, are depicted in Fig 1. The behavior 
of the order parameter si is presented in Fig 2 for different values of u and h. 

5 Discussion 

The influence of atomic vibrations and mesoscopic disorder on the properties of an insulat- 
ing double-well optical lattice are studied. Atomic vibrations result in the appearance of 
well-defined phonon degrees of freedom, provided that the interactions between atoms are 
sufficiently long-ranged. For example, local interactions, described by a delta-function po- 
tential, typical of many rarified trapped gases, do not allow for the formation of well-defined 
phonons, since their spectrum turns out to be imaginary. But long-range interactions, such 
as existing between atoms with electric or magnetic dipoles, perfectly allow for the appear- 
ance of well-defined phonon degrees of freedom. At the present time, several atomic species 
are known, enjoying such long-range dipole interactions, for instance, ^^Cr possessing large 
magnetic moments [31,32], Rydberg atoms [33] and polar molecules [34]. 

So, the consideration of the present paper concerns this type of dipolar gases loaded into 
an optical lattice. We show that atomic vibrations in a double- well optical lattice renormalize 
the strength of atomic interactions. 

Atomic vibrations can be called microscopic fiuctuations, since they correspond to the 
oscillations of separate atoms, though interacting with each other, but the motion of atoms 
being not correlated. There exists another kind of fluctuations, when atoms move in a co- 
herent manner, locally defining an order parameter that differs from that of the surrounding 
matter. The typical size of such a coherent atomic group is intermediate between the mean 
interatomic distance and the linear system size, because of which these coherent oscillations 
are termed mesoscopic. When the latter correspond to the order parameters of different 
thermodynamic phases, these oscillations can be called heterophase. Such mesoscopic fiuc- 
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tuations can occur in various statistical systems of condensed-matter type [21-27] and in 
many macromolecular systems [35]. 

Mcsoscopic fluctuations can also occur in double- well optical lattices [17]. Here, an 
ordered state corresponds to a nonzero atomic imbalance between the wells of the double 
wells. If there are no external fields, the disordered state is characterized by zero atomic 
imbalance [17]. But in the presence of external fields, disorder corresponds to an atomic 
imbalance that is not exactly zero, but smaller than the imbalance in the ordered state. We 
have generalized the previous analysis of mcsoscopic disorder in a lattice without external 
fields to the case of nonzero external fields that can be related, e.g., to the presence of a 
trapping potential. The existence of external fields essentially changes the system properties, 
leading to the occurrence of first-order phase transitions between the states with different 
order parameters, that is, between the states with different atomic imbalance. 

It is worth noting that term (16), corresponding to an external field, arises only when, in 
addition to the lattice potential, there is a nonuniform external potential, such as trapping 
potential. To show this, let us consider an ideal optical lattice with a periodic potential 

[/(r) = [/(r + a,). 

In that case, the single-atom Hamiltonian (2) enjoys Bloch waves as its eigenstates, defined 
by the eigenproblem 

Taking for atomic orbitals the well-locahzed Wannier functions [19], we have the relation 
Then for the matrix elements (5), we obtain 

k 

Thus, these matrix elements are purely diagonal with respect to the band indices m, n. 
Therefore, the effective field (16), caused by nondiagonal elements, is strictly zero, hence 
/i = 0. 

The occurrence of mcsoscopic fluctuations, causing the appearance of mcsoscopic disorder 
in a double- well optical lattice, changes the system properties in a much more dramatic way 
than microscopic fluctuations, such as atomic vibrations. 
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Figure Captions 



Fig. 1. Transition temperature Tq in units of /: (a) as a function of h for different values 
of u; (b) as a function of u for different values of h. 

Fig. 2. Order parameter si as a function of dimensionless temperature for different values 
of the disorder parameter: (a) u — 0.1; (b) u — 0.5. The lines are marked by the numbers 
corresponding to different values of the external field: (1) h — 0; (2) h — 0.2; (3) h — 0.4; 
(4) h = 0.6; (5) h = 0.8; (6) h = 1.0. 



19 




Figure 1: Transition temperature Tq in units of /: (a) as a function of h for different values 
of u] (b) as a function of u for different values of h. 
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Figure 2: Order parameter si as a function of dimensionless temperature for different values 
of the disorder parameter: (a) u = 0.1; (b) u = 0.5. The hues are marked by the numbers 
corresponding to different values of the external field: (1) h = 0; {2) h = 0.2; (3) h = 0.4; 
(4) h = 0.6; (5) h = 0.8; (6) h = 1.0. 
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